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ABSTRACT 

We study charm D - meson decays to two pseudoscalar mesons in Cabibbo 
favored mode employing SU(3)-flavor for the nonfactorizable matrix elements. 
Using D — > Kn and D s — > KK to fix the reduced matrix elements, we obtain 
a consistent fit for i] and rf emitting decays of D and D s mesons. 
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It is now fairly established that the naive factorization model does not 
explain the data on weak hadronic decays of charm mesons. On one hand 
large N c — > oo limit, which apparently was thought to be supported by D- 
meson phenomenology [1,2], has failed to explain B- meson decays, as B- meson 
data clearly demands [3] a positive value of the (^-parameter. On the other 
hand even in D-meson decays, the two body Cabibbo favored decays of D° 
and Df involving rj and rf in their final state have proven to be problematic 
for a universal choice of a\ and a 2 [4]. Annihilation terms, if used to bridge 
the discrepancy between theory and experiment, require large form factors, 
particularly for D — > K + rj/rj' and D° — > K + rj decays [4]. Further, 
factorization also fails to relate 77/7/ + n + / p + decays with semileptonic 

decays D+ — > 77/77' + e+z/ [4,5] consistently. 

Recently, there has been a growing interest in studying nonfactorizable 
terms for weak hadronic decays of charm and bottom mesons [6]. In an earlier 
work [7], we have searched for a systematics in the nonfactorizable contribu- 
tions for various decays of D° and D + mesons involving isospin 1/2 and 3/2 
final states. We observe that the nonfactorizable isospin 1/2 and 3/2 am- 
plitudes have nearly the same ratio for D — > Kir/ Kp/ K*ir/ Kai/K* p decay 
modes. In order to realize the full impact of isospin symmetry, and to relate 
-D+-decays with those of the nonstrange charm mesons, we generalize it to the 
SU(3)-flavor symmetry. 

We analyze Cabibbo favored decays of D°, D + and mesons to two 
pseudoscalar mesons. Determining the SU(3) reduced matrix elements from 
D+ -> K°tt+ and D+ -> K°K + , we obtain a consistent fit for D° -> K+n/r]/r]' 
and — > n + rj/r}' decays. 
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We start with the effective weak Hamiltonian 

H w = G F [ci(ud)(sc) + c 2 (sd)(uc)}, (1) 

where G F = ^V ud V* s and q x q 2 = <?i7 M (l — 75)^2 represents color singlet V — A 
current and the QCD coefficients at the charm mass scale are 

ci = 1.26 ±0.04, c 2 = -0.51 ±0.05. (2) 

Separating the factorizable and nonfactorizable parts, the matrix element of 
the operator (ud)(sc) in eq. (1) between initial and final states can be written 

as 

< P 1 P 2 |(Mrf)(sc)|D > = < Pi|(iid)|0 >< P 2 |(sc)|£> > 

+ < P 1 P 2 \(ud)(sc)\D > nonfac . (3) 

Using the Fierz identity 

(ud)(sc) = ±(sd)(uc) + \j:(s\ a d)(u\ a c), (4) 
c a=i 

where giA a g 2 = ?i7^(l — 75)A a g 2 represents color octet current, the nonfactor- 
izable part of the matrix element in eq.(3) can be expanded as 

< P 1 P 2 \(ud)(sc)\D > nonfac = — < P 2 \(sd)\0 >< P!\(uc)\D > 

1 8 1 

+ - < P 1 P 2 \ J2(sX a d)(u\ a c)\D > nonfac +— < P 1 P 2 \(sd)(uc)\D > nonfac . 

a=l c 

(5) 

Performing a similar treatment to the other operator (sd)(uc) in eq.(l), the 
decay amplitude becomes 

<P 1 P 2 \H W \D> = G F [a 1 <P 1 \(ud)\0><P 2 \(sc)\D> 

+a 2 < P 2 \(sd)\0 >< P!\{uc)\D > 
±c 2 (< PiP 2 \H*\D > + < PiP 2 \Hl\D >) nmfac 



+ci(< PiP2\Hl\D > + < PiP2\Hl\D >) nonfac ], (6) 

where 

Ol,2 = Ci, 2 + — , (7) 

Hi = lj:(sX a d)(uX a c), Hi = ±£(uA'd)(5A a c); 

Z a=l Z o=l 

^ = ^(5d)(fic), Hi = ±-(ud)(sc). (8) 

Thus nonfactorizable effects arise through the Hamiltonian made up of color- 
octet currents (if® and Hi ) and also of color singlet currents ( Hi and Hi 

)■ 

Matrix elements of the first and the second terms in eq. (6) can be calculated 
using the factorization scheme [1]. These are given in Table I. So long as one 
restricts to the color singlet intermediate states, remaining terms in eq. (6) are 
ignored and one usually treats a x and a 2 as input parameters in place of using 
N c = 3 in reality. It is generally believed [1,8] that the D — > Kir decays favour 
N c — > 00 limit, i.e., 

ai « 1.26, a 2 ~ -0.51. (9) 

However, it has been shown that this does not explain all the decay modes of 
charm mesons [4,5]. For instance, the observed D° — > K°rj and D° — >■ J? ^' de- 
cay widths are considerably larger than those predicted in the spectator quark 
model. Also in D — > mode, measured branching ratios for D° — > K*°r), 
— > 77/77' + p + , are higher than those predicted by the spectator quark dia- 
grams. For D+ — > 77/77' + 7r + , though factorization can account for substantial 
part of the measured branching ratios, it fails to relate them to correspond- 
ing semileptonic decays D+ — > 77/77' + e + u consistently [4,5]. In addition to 
the spectator quark diagram, factorizable W-exchange or W-annihilation dia- 
grams may contribute to the weak nonleptonic decays of D mesons. However, 
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for D — > PP decays, such contributions are helicity suppressed [1]. For D me- 
son decays, these are futher color-suppressed as these involve QCD coefficient 
c 2 , whereas for D+ — > PP decays these vanish [4] due to the conserved vector 
(CVC) nature of isovector current (ud). Therefore, it is desirable to investigate 
nonfactorizable contributions more seriously. 

It is well known that nonfactorizable terms cannot be determined unam- 
biguiously without making some assumptions [6] as these involve nonperturba- 
tive effects arising due to soft-gluon exchange. We thus employ SU(3)-flavor- 
symmetry [9] to handle these matrix elements. In the SU(3) framework, the 
weak Hamiltonians H^, H^, and for Cabibbo-enhanced mode behave 
like iPf 3 component of 6* and 15 representations of the SU(3). Since and 
transform into each other under interchange of u and s quarks, which forms 
V-spin subgroup of the SU(3), we assume the reduced amplitudes to follow 



Then, the matrix elements < PiP 2 \H^\D > can be considered as 
weak spurion + D — > P + P scattering process, whose general structure can 
be written as 



< P 1 P 2 \Hl\D > = h{P?P c m P b )Hl c] + d^P^H^ 

+e 1 (P b m P c a P m )H^ c) + h{P™P b a P c )H^ c) (11) 



where P a denotes triplet of D-mesons P a = (D°, D + , D+) and P fe a denotes 
3 ® 3 matrix of uncharmed pseudoscalar mesons, 



< P 1 P 2 \\H*\\D >=< P 1 P 2 \\Hl\\D > . 



(10) 



/ Pi 7T+ K + \ 

it- Pi K° 
\K- K° P 3 3 J 



(12) 



with 



Particle data group [10] defines the physical 77 — 7/ mixing as 

rj = T)$ cos <p — 1] sm 0; 

rf = rj 8 sin (f) + r] cos (j), (13) 

where = —10° and = —19° follow from the quadratic mass formula and 
the two photon decays widths respectively [10]. We employ the following basis 
[4] 

77 = —j=(uu + dd) sm.6 — (ss) cos 9, 
v2 

rf = —=(uu + dd) cos 9 + (ss) sin 9, (14) 
V2 

where 9 is given by 

9 = 9 ldeal -<p. (15) 
Performing a similar treatment for and H^, i.e. 

<P 1 P 2 ||^||L>> = <P 1 P 2 \\H^\\D>, (16) 

the matrix elements < PiP^H^D > are obtained from 

< P,P 2 \Hl\D > = b 2 (PTP c m P b )H^ c] + d 2 (P™P c m P b )H^ c) 

+e 2 (P b m P c a P m )H? h , c) + f2(P™P b a P c )H? b>c) (17) 

Since the C.G. coefficients appearing in the eqs. (11) and (17) are the same, 
the unknown reduced amplitudes get combined as 

b = b 1 + b 2 , d = di + d 2 , e = e x + e 2 , f = fi + /b, (18) 
when the matrix elements are substituted in eq.(6). 



There exists a straight correspondence between the terms appearing in (11) 
and (17) and various quark level processes. The first two terms, involving 
the coefficients b's and d's, represent W-annihilation or W-exchange diagrams. 
Notice that unlike factorizable W-exchange or W-annihilation diagrams, these 
diagrams are not suppressed on the basis of the helicity arguments due to 
the involvement of gluons. The third term, having coefficient e's, represents 
spectator quark like diagram where the uncharmed quark in the parent D- 
meson flows into one of the final state mesons. The last term is like a hair-pin 
diagram, where qq generated in the process hadronizes to one of the final state 
mesons. Thus obtained nonfactorizable contributions to various D — > PP 
decays are given in Table II. 

Now we proceed to determine the SU(3) reduced amplitudes b, d, e, /. First, 
we calculate the factorizable contributions to various decays using N c = 3, 
which yields 

cii = 1.09, a 2 = -0.09 (19) 
For the form factors, we use 

F DX (0) = 0.76, F^(0) = 0.83, (20) 

as guided by the semileptonic decays [8, 12], and 

F Dv (0) = 0.68, F Dr >'(0) = 0.65, 

F° sV (0) = 0.72, F° sV '{0) = 0.70, (21) 

from the BSW model [1]. Numerical values of the factorizbale amplitudes are 
given in col (iii) of Table I. 

D — > Kn decays involve elastic final state interactions (FSI) whereas the 
remaining decays are not affected by them. As a result, the isospin amplitudes 
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1/2 and 3/2 appearing in D — > i^7r decays develop different phases; 

A( J D°^K-7T + ) = -L[A 3/2 e l<5 3/ 2 + v ^A 1/2 e^], 

A(D°^K°ir°) = -±=[V2A 3/2 e iSs / 2 - A 1/2 e iS ^}, 

V3 

A{D + ^K\ + ) = V3A 3/2 e t5 ^. (22) 

which yield the following phase independent [7,11] expressions: 

\A(D°^K-7r + )\ 2 + \A(D ^K°n )\ 2 = \A 1/2 \ 2 + \A 3/2 \ 2 , 

\A(D + ^ K°7T + )\ 2 = 3\A 3/2 \ 2 . (23) 

These relations allow one to work without the phases. Writing the total decay 
amplitude as sum of factorizable and nonfactorizable parts 

A(D — > Kir) = A f (D — > Kn) + A n $ (D — > Kir), (24) 

we obtain 

A\ s l2 = -^={V2A nf (D° ^ K-7r+)-A nf (D° ^ K°ir )}, (25) 

A$ 2 = -^{A nf (D° ^ K-n + ) + V2A nf (D° ^ K\ )}, 

= ^={A^(D+^R\+)}. (26) 

The last relation (26) leads to the following constraint: 

b + d ci + Co / \ 

= 1 2 = - 0.424 ±0.042. 27 

e c 2 - ci 

Experimental value £(£>+ -> K°7r + ) = 2.74 ± 0.29% yields, up to a scale 
factor G F , 

e = - 0.094 ± 0.027 GeV 3 . (28) 
This in turn predicts sum of the branching ratios of D° — > Kir decay modes, 

B(D° ^ R-7T+) + B(D° ^ K\°) = 6.30 ±0.67% (6.06 ± 0.30% Expt.) 

(29) 



in good agreement with experiment. Using the experimental value of B(D+ — > 
K°K + ) = 3.5 ± 0.7%, we find (in GeV 3 ) 

b = +0.080 ±0.026, (30) 

d = -0.040 ±0.026. (31) 

Note that the unknown reduced amplitude / appears only in decays involving 
i] and rj' in the final state. We find that experimental values of these decay 
rates require (in GeV 3 ): 



f = 


-0.145 ±0.077 


for 


D° - 






f = 


-0.115 ±0.012 


for 


D°- 






f = 


-0.104 ±0.163 


for 


Dt- 






f = 


-0.081 ±0.073 


for 




-> V 'ti + . 


(32) 



In Tables III, we calculate branching ratios for all the four rj, emitting decay 
modes for different choice of /, for = — 10° and —19°. It is clear that for / = 
—0.12 and <fi = —10°, all the branching ratios match well with experiment. For 
the sake of comparison with factorizable terms, nonfactorizable contributions 
to various modes for / = — 0.12 are given in column (iii) of the Table II. 
Color-suppressed decays obviously require large nonfactorizable contributions. 
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Table I 



Spectator-quark decay amplitudes ( x G f GeV 3 ) 



Process 



D+ -»• K y h 



Amplitude 



10 L 



19 L 



D°- 
D°- 

D° 

D° - 

Dt- 

Dt 

Dt 



* K°T] 
> 7T°7r + 

-> ?77r + 



+ a 2 f K (ml-ml)F ^(mj c ) 



aifn(m 2 D - m 2 K )F ( 



DK i 



?Dtt, 



-^a 2 f K (m 2 D -m 2 w )F^ 
-^a 2 sin9f K (m 2 D - m^F^ 
-^a 2 cos9f K 



rn 



(m 



D 



rn 2 ,)F^\m 2 K ) 



a 2 f K (m 2 Da 



m 




)F^ K {m 2 K ) 



K 



-a 1 cos6f n (m 2 Dg 

,2 ^2^%'^ 



aisinO U(m 2 Ds - m;,)F 



+0.311 

+0.354 
-0.030 

-0.016 

-0.013 

-0.035 


-0.261 
+0.213 



+0.311 

+0.354 
-0.030 

-0.019 

-0.010 

-0.035 


-0.216 
+0.243 



Table II 



Nonfactorizable contributions to D — > PP decays ( x Gf GeV 3 ) 



Process 


Amplitude 


0=- -10 u 


= -19 u 


D + ^ 


2(ci + c 2 ) e 


-0.141 


-0.141 


D° -> K-TT+ 
D° -> K°n° 
D° -> i^ ?7 
D° -> 


c 2 (6 + d + e) 
(-6 - d + e) 
Cl [^ (6 + rf + e + 2/) cos#(6 + d + /)] 
Cl [^ (b + d + e + 2f) + sm0(fc + d + /)] 


+0.028 
-0.119 
-0.115 
-0.256 


+0.028 
-0.119 
-0.154 
-0.235 


Z>+ -> 7T 7T + 
£> + -> ^7T + 
£>+ -> t/tT + 


ci (-6 + d + e) 


c 2 [V2sin6 (-b + d + f) - cos6(e + f)} 
c 2 [V2cos9 (-6 + d + /) + sira0(e + /)j 


-0.268 


+0.046 
+0.199 


-0.268 


+0.076 
+0.189 
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Table III 



Branching 


:% 


i of r]/r]' emitting decays 


includin 


g nonfactorization terms 


Decay 






> = -10° 








= -19° 




Expt. 




/ 


= -0.10, 


-0.12, 


-0.14 


/ 


= -0.10, 


-0.12 


-0.14 




D° -> r)K° 




0.53 


0.59 


0.66 




0.86 


1.02 


1.19 


0.68±0.11 


D° -» rfK 




1.28 


1.81 


2.43 




1.04 


1.51 


2.06 


1.66±0.29 


D+ -» »)7r+ 




1.93 


1.87 


1.82 




0.86 


0.80 


0.73 


1.9±0.4 






5.17 


5.64 


6.13 




5.73 


6.22 


6.72 


4.7±1.4 
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